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SUMMARY 

The alternating direction implicit (ADI) method is adopted, modified, and applied to the Reynolds 
equation for thin, gas fluid films. An efficient code is developed to predict both the steady-state and 
dynamic performance of an aerodynamic journal bearing. An alternative approach is shown for hybrid 
journal gas bearings by using Liebmann’s iterative solution (LIS) for elliptic, partial differential 
equations. 

The results are compared with known design criteria from experimental data. The developed 
methods show good accuracy and very short computer running time in comparison with methods based 
on an inverting of a matrix. The computer codes need a small amount of memory and can be run on 
either personal computers or on mainframe systems. 


INTRODUCTION 

The Reynolds partial differential pressure equation for a compressible, thin, fluid film is nonlinear 
and involves a great deal of numerical effort to solve. From the early sixties to the present time several 
methods were developed for numerically solving this equation. Major contributions were made by 
Castelli and Pirvics, Booy, Coleman, Castelli in Gross’ book, and others (refs. 1 to 4). The columnwise 
influence coefficients method introduced by Castelli was one of the best at that time. The present work is 
focused basically on two methods (alternating direction implicit method (ADI) and Liebmann’s iterative 
solution (LIS)), which add important advantages to the numerical solution technique of the compressible 
Reynolds equation. 

The Reynolds equation in dimensionless form for a compressible-fluid-film journal bearing subject 
to isothermal conditions is (ref. 5) 

1 (h 3 £e!) + 1 (h 3 *!) = 2A + 2A, + Ufi ^ (1) 

d0 do di dz dO dz dr 

Both steady-state and dynamic analyses to calculate the static and dynamic performance of an 
aerodynamic journal bearing are shown in the appendix A. Starting from equation (1), the following 
partial differential equation are developed in appendix A: 
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Equation (2) is used to perform the steady-state characteristics; while equations (3) and (4) are used to 
compute the dynamic performance. The numerical integration techniques for these equations are 
developed and discussed in herein. 


SYMBOLS 


a ij’ kjj> c ij> 


djji e ij> r jj 


B c 

B ij 
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"D 
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coefficients of grid node eq. (6) 
supply orifice radius, m 

symbolic damping coefficient used in stability analysis (see eq. (A-27)) 

dynamic damping coefficient, N-m _1 -s -1 

journal bearing radial clearance, m 

orifice supply discharge coefficient 

journal bearing diameter, m 

orifice supply hole diameter, m 
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K ij 

K c , K c0 
L . 
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M c 

M,N 


eccentricity (see fig. l), m 

maximum iteration step error (see eq. (12)) 

load capacity (the resulting force of the pressure distribution), N 

F/(p a LD) dimensionless load capacity (eq. (18)) 

component of F along and perpendicular to the centerline, respectively (see fig. 1) 
steady-state component of F 

dynamic component of F along x and y directions (see fig. 1) 
whirl frequency ratio, i//ft 
unstable whirl frequency ratio, i/ 0 /fi 

dimensionless supply flow part that is function of p m /p g (see eq. (B-ll)) 

dimensionless film thickness, TT/C 

film thickness, m 

steady-state component of h 

film thickness at the supply orifice (see fig. 13) 

— I 1 / 2 the imaginary unit 

the unit vector along Xj direction (see fig. 14) 

stiffness, N*m — 1 

dimensionless stiffness (see eq. (21)) 
dynamic stiffness coefficient 

symbolic stiffness coefficient used in stability analysis (see eqs. (A-27) and (A-29)) 
bearing length, m 

length vector part of the C m boundary of integral (B-7) (see also fig. 14) 
rotor mass allocated to one bearing; for a symmetric rotor M is half of the rotor mass, kg 
corresponding rotor mass, allocated to one bearing, required to make the bearing unstable, kg 
number of grid points in i and j directions 
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n 


P 

P 

Pa 

Pl»’ ?2a 
Pm 
Po 
Ps 

Pi’ P2 
Q 


Qo’ Qi> Qj 

R 

P 

r, t 
S T 
T 
t 

V, 

w 


x, y 



z 


unit vector outward, perpendicular to the C m boundary (see n- in fig. 14) 

total number of supply orifices 

pressure, Pa 

dimensionless pressure 

ambient pressure, Pa 

boundary pressure at bearing edges 

pressure downstream of the supply system 

steady-state component of the pressure p (see eq. (A-6)) 

supply pressure, Pa 

perturbation components of pressure p (see eq. (A-6)) 
dimensionless pressure variable, (ph) 2 

steady-state and perturbation components of Q (see eq. (A-7)) 

dimensionless flow through the pocketed orifice supply restrictor system 

journal bearing radius, m 

gas constant (see eq. (B-9)), J*kg **K 1 

coordinates (fig. 1) 

threshold speed, ft(M c C/W) 1 / 2 

absolute temperature (see eq. (B-9)), K 

time, s 

axial velocity (shaft speed along z direction), m*s 1 
bearing load, N 

whirl amplitude of the journal center in bearing (see eqs. (A-20) and (A-25)) 
fluid film coordinates (see eq. (A-l) and fig. 14) 
impedance for translatory motion 
axial coordinate parallel to rotor axis 
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7 

adiabatic exponent 

Az 

increment in z direction 

a e 

increment in 0 direction 

6 

pocketed orifice compensation factor, a 2 /dC 

e 

eccentricity ratio, e/R 

( 0 

eccentricity ratio under static load 

( 0^1 

dimensionless tangential whirl amplitude 

( 1 

dimensionless radial whirl amplitude 

9 

angular coordinate from centerline (see fig. 1) 

A 

bearing number (see eq. (A-2) and (B-5)) 

K 

equivalent bearing number in z direction (see eq. (A-2)) 

K 

restrictor coefficient (see eq. (B-9)) 

A 

relaxation coefficient (see eq. (10)) 

9 

o 

dynamic viscosity, N*s-m 

V 

whirl frequency, rad-s~ 1 

"o 

unstable whirl frequency 

T 

dimensionless time (see eq. (A-l)) 

* 

attitude angle, deg 

^0 

attitude angle under static load, deg 

n 

rotation frequency, rad-s^ 1 


shaft rotation frequency, rad-s” 1 

n 2 

bearing housing rotation frequency, rad-s _1 

Subscripts: 


ij 

grid index 

r 

centerline direction (see fig. 1) 
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t perpendicular to center-line direction (see fig. 1) 

x x-direction (direction of the static load on bearing, see fig. 1) 

y y-direction perpendicular to x (see fig. l) 

Superscripts: 

n old iteration step 

n + 1 new iteration step 

(n + l)/2 half step iteration between n and n + 1 


NUMERICAL COMPUTATION 

In order to integrate equations (2) to (4), a finite-difference numerical method is adopted. This 
method can handle a plain, cylindrical, journal bearing of finite length (fig. 1) with a high degree of 
accuracy. Using a rectangular grid with constant spacing between nodes in both directions, the partial 
derivatives become 
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At node (ij) of the grid, equation (2) has the general form 


+ bi j Q °i+u + C 'j Qo ij-i + di J^°ij+i + ei J^°ij ri J 


( 10 ) 


Thus, a linear system of M x N equations results if the grid has M and N points in the z and 6 direc- 
tions, respectively. Unfortunately, the coefficients of equation (10) depend on the value of the variable Q 
so that an iterative method has to be used to compute the steady-state component of Q. 

First, the system of equations (10) must be solved. Then, since 

Qo - Po (“) 


and all the steady-state characteristics, the equations (3) and (4) can be solved, using the same finite 
difference method, to calculate the dynamic performance. 


ALTERNATING DIRECTION IMPLICIT SCHEME 

The linear system (10) can be efficiently solved using the alternating direction implicit (ADI) scheme 
(ref. 6), which changes the linear system of equation (10) into a tri-diagonal linear system using a proper 
iterative manner. The tri-diagonal system can be solved with Thomas’ algorithm. In order to use the 
ADI method, equation (10) was written as follows for the first and second half step of the iteration n + 1: 
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This same scheme was applied to equations (3) and (4) to obtain the dynamic characteristics. 

Using the ADI method, the computer running time for one case (one eccentricity ratio) is reduced 
from 1 hr to 1 to 5 min (using a 20x9 grid and a 386, 20-MHz personal computer). To decrease the total 
time to reach solution convergence, equation (11) can be used alone. 

The convergence is difficult for equation (2), (3), or (4) for bearing numbers (eq. (A-2)) greater than 
10. Also, the convergence of the numerical solution is sometimes difficult when the stability of the 
bearing is analyzed, especially at high bearing numbers. One possibility to improve the convergence is to 
use a relaxation technique: 


Q n+1 - 


AQ n+1 + ( 1 - A) Q n 


(14) 



where 


A = 0 — 1 

A second possibility is to convert the dimensions of the grid from a large one (e.g., from 
NxM = 40 x 33 or 30x25 or 20x17), used to compute the steady-state component Q 0 of the variable Q 

(eq. (A-7)) and the derivatives of Q 0 and h 0 , to a small one (e.g., 10x9) and then to compute the 
dynamic components Qj and Q 2 of the variable Q. 


LIEBMANN’S ITERATIVE SOLUTION 

The discretization of the Reynolds equation results in a large linear system of equations that needs 
to be solved. For example, a 20 by 9 grid involves 180 linear algebraic equations like equation (10). 
There are a maximum of five unknown terms per line in equation (10). For a large-sized grid, many of 
the terms will be zero. For such a sparse system, the elimination methods waste significant computer 
memory to store zeros. To avoid these difficulties, the LIS method (ref. 7) can be used as an alternative 
to ADI. In this technique equation (10) is expressed as 

_ R ij “ A ij Qi-lj “ B ij Qi+lj ~ C ij Qij-1 ~ D ij Qij+1 (J5 


and solved iteratively for j = 1 to N and i = 1 to M. Because equation (10) is diagonally dominant 
for most of the bearing running conditions, this procedure will eventually converge to a stable solution. 
Note that the coefficients C, D, and E depend on the value of the variable Q from the previous 
iteration. As with the ADI method, a relaxation technique can be used to improve the convergence of 
this method. The relaxation coefficient, A, can be variable and under control of the maximum error for 
each iteration step calculated as 


w = Max 


(%)" +1 - (Qij)" 


(Qij)" +1 


(16) 


for 


i = 1 — ► M , j = 1 — N 


(17) 


One of the most difficult cases regarding the convergence is the hybrid, compressible-fluid, film 
bearing. This case balances the flow through the supply restriction system with the flow through the 
fluid film to calculate the downstream pressure in the supply port. This requires an additional iteration 
loop. Also, the equation of the flow through a restriction Bystem, like a pocketed orifice, is a combination 
of four relations that is a function of the type of flow (choked or unchoked) and the direction of the flow 
(in or out). The analysis of a hybrid, compressible fluid bearing is outlined in appendix B. 
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RESULTS AND DISCUSSION 


ADI Method 

Both steady-state and dynamic performance results over a range of bearing number A (A-3) from 
0.903 to 9.03 and an eccentricity ratio between 0.1 and 0.6 are very close to known computed data from 
(refs. 8 and 9). Figures 2 and 3 show the basic steady-state characteristics: load capacity and attitude 
angle (fig. 1 and eq. (A-19)) versus bearing number. The dimensionless form of the load capacity is 

F = - - (18) 

P a LD 


The numerically computed steady-state data using the ADI method agree well with Raimondi data 
(ref. 8) for both grid dimensions used (20x9 and 10x5). The same conclusion arises in figures 4 and 5 
where the whirl frequency ratio and threshold speed criteria are shown. The whirl frequency ratio is 


f 0 



(19) 


and the threshold speed is 
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The notation used in equations (19) and (20) are described in appendix A. In figures 4 and 5 the ADI 
computed data are compared with Chang and Pan’s computed data (ref. 9). An exception to this general 
conclusion may be for the threshold speed at low bearing numbers and especially at small eccentricity 
ratios (left bottom corner of fig. 5) where the ADI computed data show more sensitive behavior to the 
“half frequency instability” than the data from Chang and Pan. No experimental data were found for 
this particular range of input data. 


LIS Method 

The LIS method is first used to compute the steady-state performance of an aerodynamic journal 
bearing. To assess the accuracy of the LIS method, both LIS and Raimondi data are shown in figure 6 
for load capacity and in figure 7 for attitude angle. Very good agreement was obtained for a range of the 
bearing numbers between 0.1 and 17. The grid has 30x9, 40x9, and 60x9 nodes, and the computing 
time for one case (one point in either fig. 6 or 7) is regularly less than 1 min. 

The LIS hybrid gas bearing code was first used to compute the steady-state characteristics of a 
double-plane supply aerostatic journal bearing with a length to diameter ratio L/D of 2. The output 
data for the dimensionless stiffness as was defined in (ref. 10) as 


9 



( 21 ) 


K = 1 + 6 1 C K 

1 + 2 S 2 (p g - Pa) LD 
3 

versus restrictor coefficient (eq. (B-9)) was plotted in figure 8. The notation used in equation (21) is 
shown in appendix B. Based on an average stiffness calculation, the eccentricity of the shaft was 0.05, 

0.1, and 0.5 of the radial clearance (curves LIS 0.05, 0.1, and 0.5 respectively). For comparison, the 
corresponding curve from Wilcock (ref. 10) and Pink’s experimental data (ref. 11) are also plotted in 
figure 8. Pink’s data (ref. 11) were obtained using the stiffness near zero eccentricity and the load at 0.5 
eccentricity. The LIS computed data agree well with Pink’s experimental data. 

The computed total load and stiffness for a hybrid journal gas bearing is shown in figures 9 and 10. 
Both externally pressurized and dynamic effects contribute to pressure rises in the fluid film. The 
dynamic effects contribute a significant amount to the total load, especially at high eccentricity ratios 
and high speeds. Figure 10 and especially figure 11 indicate that, for small eccentricity ratios (near zero), 
the dynamic effects can diminish the bearing performance in turning the shaft at low speeds. The radial 
flow due to the rotation of the shaft appears to change the externally supplied flow. This change results 
in the load and stiffness of the hybrid bearing (turning speeds of 5000 and 15000 rpm) being lower than 
the aerostatic bearing curve (0 rpm) for eccentricity ratios 0 to 0.2. This phenomenon is controlled by 
the turning speed and by the radial clearance ratio. Figure 11 shows how radial clearance affects the 
bearing stiffness for turning speeds of 0 and 5000 rpm. 


CONCLUSIONS 

Two methods (alternating direction implicit (ADI) and Ziebmann’s iterative solution (LIS)) were 
adopted and modified to integrate the compressible pressure-differential equation. The ADI method was 
used to compute the steady-state and the dynamic performance of an aerodynamic journal bearing. Both 
steady-state and dynamic results are in good agreement with known design criteria. The LIS method 
was used to calculate the pressure distribution of a hybrid journal gas bearing. Because a hybrid bearing 
needs an additional iteration loop to balance the flow through the supply system with the flow through 
the lubrication film, convergence takes longer. Both numerical methods (ADI and LIS) are more 
efficient, in terms of running time of the codes, than other known methods. However, the ADI method 
sometime shows floor convergence, especially for dynamic computation at very low or very high bearing 
numbers. The LIS method can handle the aerodynamic, the aerostatic, and the hybrid bearings and can 
generate useful data, especially for hybrid gas bearings. Furthermore, the LIS method can provide data 
with minimal computer time. Both methods (ADI and LIS) provide data with running computing times 
between 10 to 50 times faster than other methods that involve inverting a matrix. The running time for 
the ADI or LIS method is proportional to the grid dimension. 
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APPENDIX A. — STEADY-STATE AND DYNAMIC ANALYSES OF A PLAIN, CYLINDRICAL, 
FINITE LENGTH, SELF-ACTING GAS JOURNAL BEARING. 


The Reynolds equation in dimensionless form for a compressible-fluid-film journal bearing subject to 
isothermal conditions is (refs. 4 and 5) 


where 
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f = l (n = n, + n 2 ) 


(A-3) 


and where flj is the shaft rotating speed, ft 2 is the bearing rotating speed, v is the whirl frequency (in 
radians per second), v z is the shaft axial speed, and /x is the dynamic viscosity. 

Both circumferential and axial movement are taken into account in equation (A-l), and a complex 
form for time-dependent terms is used for the dynamic analyses. 

In order to integrate equation (1), a natural linearization procedure is applied (ref. 1) by writing 
Reynolds equation in terms of a new variable Q = (ph) 2 . For the dynamic analyses, a small perturba- 
tion technique was used. Here, the journal center was assigned a static eccentricity ratio e 0 and a 
corresponding attitude angle <f> G . This position was perturbed by a small amplitude harmonic motion: 


e = f 0 + < 1 e r , <p = <p Q -I- p 1 e r 


(A-4) 


The corresponding perturbed forms of the dimensionless film thickness and pressure can be written as 


h = h 0 + €je T cos 9 + sin 9 


(A-5) 


P = P 0 + f l ef Pi + P2 


(A-6) 
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Introducing a new variable 


Q = (ph) 2 = Qo + 2 «i eT Qi + 2 Wl*' Q 2 

and substituting equations (A-5), (A-6), and (A-7) into equation (1) yield 
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The equation (2) is the conventional steady-state equation. It is nonlinear and was solved in finite- 
difference form by an iterative method. The last two equations depend on the steady-state solution. 
They are linear equations with complex dependent variables and were solved numerically. 

Taking into account equations (A-5), (A-6), and (A-7) the following relations can be established: 


Qo ” Po 2 N) 2 

Qi - Po 2h o cos 9 + Po h o 2 Pi 
Q 2 - Po 2h O sin 9 + Po h 0 2 P 2 
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or 
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_Qi_ 
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ho 
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respectively. The boundary conditions are 


2 2 

Po = Pla> Qo = Pia^o 
z = 0, Pi =0, Qj = pj a h 0 cos0 

P 2 = °> Q 2 = Pu h 0 sin * 


2 2 
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L 2 

z = — > Pi = 0, Qi = p 2a h 0 cos* 
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P 2 =0, Q 2 = P 2a h 0 sin* 


For a journal bearing a periodic condition is also necessary: 


Pfl=0 - Pfl=2r 


(A-9) 


(A- 10) 


(A- 11) 


The resulting hydrodynamic force will be computed by integration of the film pressure. The force will 
have a radial component in the opposite direction of the eccentricity (on the centerline) and a tangential 
component in the direction of the attitude angle under static load: 


F r 
Ft 

A first-order Taylor series expansion of the force components around the steady-state equilibrium position 
can be written in the form 


f = P a R2 J J (Po - 1 + Pl«l eT + P 2 f oPl eT ){ 8 i° S /f 
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(A- 13) 


where 


F r = F rO + Z rr f l eT + Z rt e 0 P\ eT 
F t = F tO + Z tr e l eT + Z tt € 0 'Px ^ 

Z ij = Kjj + i v By 
i = r,t; j = r,t 


(A-14) 


are the bearing impedances. From equations (A-12) and (A-13) the steady-state force components and 
the impedance coefficients are 
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The steady-state force is 


F = fiio 1 + F, 0 J ; W = -F 


(A-18) 


where F is the load capacity of the bearing and W is the bearing load. The angle between the force F 
and the centerline axis is defined by the equation 


* 



(A- 19) 


Under dynamic conditions, the journal (shaft) center whirls in an orbit around its static equilibrium 
position. The corresponding dynamic bearing reaction is actually a nonlinear function of the amplitude 
and depends implicitly on time. In a thorough analysis it is necessary to consider the rotor and the 
bearing simultaneously as was done, for example, in the reference 12. In most practical situations, the 
amplitude of the shaft motion is, of necessity, rather small. In these cases a linearization of the bearing 
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reaction is permissible (ref. 5). Then it becomes possible to treat the bearing separately and to represent 
the bearing reaction by means of coefficients: 


K = - - K„y - 

By = - Kyx x - - K,,y - _ 


(A-20) 


Equation (A-20) is only valid when the journal motion is harmonic, and when 


x = x e tvt = x e r 
y = y e tvt = y e T 


(A-21) 


The equation (A-20) can be written in complex form as 


F x = - V 


F = 

y 


Z yx x 


Z xy^ 

Zyyy 


(A-22) 


where 


^ij ” ^ij + i ^ By 
i = x,y; j = x,y 


(A-23) 


For a given bearing geometry the dynamic coefficients apply to a given static load on the bearing and a 
given rotor speed. The dynamic coefficients also depend on the whirl frequency, and they are actually 
impedances of the gas film. 

Note that the x-axis was chosen along the direction of the steady-state load. In this case, the 
dynamic coefficients related to the x-y axis (fig. 1) can be calculated with the equations 

C = (Kjj. cos<£ + K tr sin^)cos^ - (K^ cos^ + K tt sin^)sin^ 

C K = (Kjj. cos^ + K tr sin^)sin^ + (K rt cos^ + K tt sin^)cos^ 

(A-24) 

C = (K^ sin^ — K tr cos^)cos^ — (K^ sin^ — K tt cos<£)sin^ 

C Kyy = (K^ sin<f> - K tr cos^)sin^ + (K rt sin^ - K tt cos^)cos^ 

Three parameters are required to specify the coefficients. The most common rotor calculations, 
other than those for stability investigations, are for critical speeds and unbalanced response. For those 
conditions the journal whirls synchronously with the rotation such that the whirl frequency equals the 
rotating frequency. 
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“ = (B^ cos (ft + B tr sin^)cos^ - (B rt cos^ + B tt sin^)sin^ 

i/" 


= (B^ cos<j> + B tr sin^)sin^> + (B rt cos^ + B tt sin^)cos^ 

l/u 


C B yx 

~w 


(Bjj sin^ - B tr cos(j>)cos<l> — 


(B rt sin^ - B tt cos^)sin^ 


(A-50) 


JL = (Bp,, sin^ - B tr cos^)sin^ + (B rt sin^ - B tt cos <f> )cos <f> 

l/u 

In a stability calculation, it is necessary to evaluate the bearing coefficients over a frequency range, 
usually around one half of the rotating frequency. On the basis of these values, a stability analysis can be 
performed. If M is the rotor mass supported by the bearing (M = 1/2 of the rotor mass when the rotor is 
symmetric and rigid) and the bearing is represented by its four Z impedance coefficients, the motion 
equation can be written as 


(Z„ - M S) 


xy 


Zy X (z yy - M v l ) 


= 0 


(A-26) 


The threshold of instability occurs when the determinant of the matrix is zero. Noting that 

Z = K c + ii/B c , K c = M if 2 


(A-27) 


we can solve the determinant equation to get 


Z = i (Z„ + V - 


N 4 


- <Z= - v>’ + Z .T z 


xy yx 


(A-28) 


(For stability calculations only the solution with a negative sign in front of the square root proves to be 
of interest (ref. 5).) At the threshold of instability, Z must be real. The imaginary part of Z can be 
evaluated over a range of frequencies to find that frequency value v Q at which the imaginary part of Z is 
zero. The corresponding mass is the mass required to make the bearing unstable under the selected 
working conditions and is 

M,. = — (A-29) 


16 



APPENDIX B— HYBRID GAS JOURNAL BEARING ANALYSIS 


The geometry and notation for a double-plane, externally pressurized journal bearing is shown in 
figure 12. The Reynolds equation for steady-state condition of a gas journal bearing becomes 


a 

h 3 dp 
—P— 

+ R 3 -^ 

h 3 ap 
—p— 

= 6R 2 n 

(B-l) 

d6 

P dO J 

di 


ee 



Supposing an isothermal evolution of the gas such that 


- = constant 


(B-2) 


and using the new natural linearization variable 


Q = (Ph) 2 


(B-3) 


equation (B-l) may be written as 


aV n+1 ) ^Q^ 1 ) 


ar 


dz* 


iah A 1 

ha# h 


V'qWJ 


aQ( n+1 ) 


d$ 


2 a 2 h A 1 aQ^ 


h d $ 2 2 h (q(n)f /2 90 J 


q( h+1 ) 


a i aQ^ n ) 


2b 




e» 


(B-4) 


where the bearing number is 


A = 


XI 

R 

A 


(B-5) 


Equation (B-4) is written for two successive iterations. To improve the convergence due to the nonlinear 
terms of the compressible Reynolds equation, a Newton-Raphson method 


f( n +l) 

0*9 

= f( n ) 

Q,— 
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farf 


axj 


axj 
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AQ + 


/ \ 
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af 

(n) A 

aQ 


f \ 


3x 

d 

aQ 


V. / 

< 





is used to write equation (B-4). 


(B-6) 
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Equation (B-4) is valid for all points covering a bearing grid except the edges and the fluid supply 
points. For a pocketed orifice restrictor supply system (fig. 13), the flow through the supply restrictor 
and through the film has to be balanced to calculate the pressure downstream of the supply restrictor: 

j ( Aph - Rh 3 p Vp ) n dl = q m (B-7) 


The left side of the equation (B-7) represents the flow through the film. For computational convenience 
the supply hole geometry is slightly modified as shown in figure 14. The right side of the equation (B-7) 
is the flow through the supply orifice restrictor: 




2*-A g p 8 


1 + 6 2 

h 2 + S 2 

m 




m 


(B-8) 


where the restrictor coefficient is 


and 


A,= 


6/iiiQa 2 


PsC" N 


5RT 

1 + S 2 


(B-9) 


6 = 


d C 


(B-10) 


G m - C D 



1 


7 2 

• 2 7 

Pm 

7 - 1 

[Pb 


1 - 


/ XT* 

Pm 


vPh, 


' > 

7+1 

' > 

7 

2 

2(7-1) if Pm < 

2 

7-1 

U+!j 

P. 

.7+1. 



if 


’ ' 

7 

2 

^7+1, 

7 + 1 


Ps 




f \ 

2 

/ ^ 

7~1 


* 

27 

P. 

7 

i - - 

7 

>■ 


7 - 1 
> 



IPmJ 




> 1 


P. 


(B-ll) 
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(B- 12) 


The discharge coefficient C D (B-12) was established taking into account data from references 4 and 10. 

Note that p m is the pressure downstream of the supply system and p 8 is the supply pressure. 
Having a first-guess pressure distribution, the pressure p m can be calculated with equation (B-7) and 
used like a boundary condition for equation (B-4) in the next iteration. 
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Figure 3. — Aerodynamic journal bearing attitude angle versus 
bearing number. ADI method results. 



Eccentricity 

ratio 

Grid size 

Source 

* 

0.1 

10x5 ^ 

J 

□ 

.4 

10x5 

> ADI 

1 

A 

.1 

20x9 

0 

.4 

20x9 J 


— 

.1 

— 1 

L Raimondi 


.4 

— j 

' (ref. 8) 




Eccentricity 

ratio 

Grid size 

Source 

* 

0.1 

10x5 "N 


□ 

.4 

10x5 l 

ADi 

A 

.1 

20x9 f 

0 

.4 

20x9 J 


.1 


Chang and 

-+■ 

.4 

—I 

Pan (ref. 9) 



Figure 2. — Aerodynamic journal bearing dimensionless load Figure 4. — Aerodynamic journal bearing whirl frequency 

versus bearing number ADI method results. ratio versus bearing number. ADI method results. 
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Figure 5,— Aerodynamic journal bearing threshold speed Figure 7.— Aerodynamic joumaJ bearing attitude angle 

versus bearing number. ADI method results. versus bearing number. US data. 
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Figure 6.— Aerodynamic journal bearing dimensionless 
load versus bearing number. US data. 
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Figure 8.— Comparison of US computed stiffness with design 
data and experimental data. Aerostatic journal bearing; 
pressure ratio Pg/P„ 5; length to diameter ratio, UD, 2; 
double supply plane; zero speed. 
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Figure 9. — Total load versus eccentric rty ratio as function of 
turning speed. Hybrid air bearing; length, L, 100 mm; 
diameter, D, 50 mm; radial clearance, C, 0.02 mm; double 
supply plane; supply pressure, p*. 0.5 MPa. 
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Figure 1 0. — Stiffness versus eccentricity ratio as function of 
turning speeds. Hybrid air bearing; length, L, 100 mm; 
diameter, D, 50 mm; radial clearance, C, 0.02 mm; double 
supply plane; supply pressure, p*, 0.5 MPa. 



Radial 

Turning 

clearance, 



mm 

rpm 


0.020 

0 

-0- 

.020 

5000 


.015 

0 

♦- 

.015 

5000 

— At— 

.010 

0 


.010 

5000 



Figure 11. — Stiffness versus eccentricity ratio as function 
of radial clearance at two turning speeds. Hybrid journal 
bearing; length, L, 100 mm; diameter, D, 50 mm; double 
supply plane; supply pressure, p* 0.05 MPa 
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